CLNS 98/1551 
CLEO 98-6 



A limit on the mass of the v T 



CLEO Collaboration 
(February 7, 2008) 



Abstract 

A limit on the mass of the tau neutrino ra Vr is derived from 4.5 x 10 6 tau 
pairs produced in an integrated luminosity of 5.0 fb^ 1 of e + e~ — > 7* — > 
t + t~ reactions at center of mass energies ~ 10.6 GeV. The measurement 
technique involves a two-dimensional extended likelihood analysis, including 
the dependence of the end-point population on m UT , and allows for the first 
time an explicit background contribution. We use the the decays r — > 5tti/ t 
and r — > 37r27r ^ T to obtain an upper limit of 30 MeV/c 2 at 95% C.L. 
PACS: 14.60Lm, 12.15Ff, 13.35Dx 
Keywords: Tau Neutrino Mass 
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Although the tau neutrino has never been directly observed, the question of its mass is an 
important issue in particle physics and cosmology. While the requirement that the density of 
primordial relic neutrinos from the Big Bang not over-close the Universe restricts the mass of a 
stable neutrino to be less than ~ 100 eV/c 2 |l|-[|, unstable neutrinos are less restricted. Big Bang 
nucleosynthesis models allow for a massive neutrino in the 10 to 31 MeV/c 2 range for lifetimes 
in the 0.01 to 40 second interval [||[]. This astrophysically allowed range currently overlaps the 
experimentally accessible bounds. The most stringent accelerator-based limits on m UT are derived 
from one-dimensional mass fits and two-dimensional energy versus mass fits of the hadrons in tau 
decays. The sensitivity of the method is largest when the invariant hadronic mass is large, and 
there is little phase space for the (unseen) neutrino such as in the decays modes r — > 57r(7r°)&v 
and t — ► 37r27r°^ r . In addition, due to its large branching ratio, the decay r — ► "ii:v T can be used 
in spite of its lower inherent sensitivity per event resulting from its low average hadronic decay 
product mass. 

The lowest published upper limitQ, by the ALEPH collaboration ||, is 24 MeV/c 2 and is ob- 
tained with a two-dimensional likelihood fit to a region near the endpoint of 25 57r(-7r ) candidate 
events. By doubling this dataset and adding in 3tt decay candidates, the ALEPH collaboration Q] 
has recently reported an upper limit of 18.2 MeV/c 2 . The DELPHI collaboration |8| has claimed 
a limit of 33 MeV/c 2 in the 37r mode from a mass fit, but notes that this limit is sensitive to the 
contribution of a possible higher mass resonance in the 3tt mass distribution. DELPHI estimates 
that the model dependence in the 37r mode decreases the limit to 62 MeV/c 2 . The OPAL col- 
laboration Q quotes a 35.3 MeV/c 2 limit from a two-dimensional fit to the missing momentum 
and missing mass in 3n decays using the event thrust axis as an estimator of the tau direction. 
Combining this with an earlier limit of 74 MeV/c 2 from 5-7T decays |10], a 29.9 MeV/c 2 limit is 



obtained. ARGUS [11,12| has published a 31 MeV/c 2 limit based on a one-dimensional mass fit to 
20 5tt events. The previous CLEO limit (b| of 32.6 MeV/c 2 was set with a fit to the mass spectrum 
of 60 5ir events, and 53 37r27r° events. 



The current analysis uses data collected with the CLEO II detector [14] at the Cornell e + e~ 
storage ring (CESR) at center of mass energies near 10.6 GeV. The data set consists of 5.0 fb^ 1 of 
integrated luminosity corresponding to 4.5 million produced tau pairs. The tau's, created mainly 
via e + e~ — > 7* — > t + t~, have the full beam energy of 5.29 GeV and recoil back to back, modulo 
initial state radiation effects. The decay modes studied are r — > §-kv t and r — > 3ir2ir u T recoiling 
against a leptonically-decaying tau (the "tag"). An upper limit on m v is obtained from a two- 
dimensional likelihood fit to hadronic mass and energy in these decays. The likelihood includes 
for the first time a background contribution, as well as a factor describing the expected number of 
events as a function of neutrino mass, and measured event errors for hadronic mass and energy. 

In order to avoid a spuriously low limit from a background event, the event selection criteria are 
chosen to obtain a very pure event sample of r — ► 5ttv t or 3tx2tx®v t , and roughly follow those of our 



last publication [13]. A clean sample of tau's is obtained by selecting events in which the tag decays 



6 Recently || it has been pointed out that these arguments do depend upon which set of astro- 
physical observations are used. 

7 A11 neutrino mass upper limits quoted herein are at the 95% confidence level. 
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leptonically: r — > ez/z/f or r — > \ivv T . We thus require the event to contain an isolated lepton (e 
or \i) in the barrel (| cos#| < 0.71) recoiling against 5 or 3 charged tracks, with event hemispheres 
defined by the event thrust axis. Exactly 6 or 4 charged tracks consistent with production at the 
nominal interaction point are required, after pairs of tracks consistent with photon conversions 
and Dalitz tt° or Kg decays are rejected. Backgrounds from two photon events are minimized 
by requiring that the total missing momentum point away from the beamline, and that the total 
energy transverse to the beam axis be above 200 MeV. To supress hadronic events, the total visible 
energy is required to be below 90% of the total center of mass energy. All hadronic tracks must 
also have a specific ionization (dE/dx) consistent with that of a pion. We correct for dE/dx energy 
loss when computing energies and invariant masses. 

Isolated showers in the calorimeter with photon-like lateral shower profiles are identified as 
photons. In the barrel, photon candidates must have an energy above 60 MeV, while in the range 
0.71 < |cos#| < 0.95, they must have an energy above 100 MeV. Exactly or 2 neutral pions 
must be reconstructed on the signal side, with none on the tag side. These 7r° candidates must 
have an opening angle of less than 60° and a mass constrained fit of the ir° candidates must have 
a confidence level greater than 5%. If the photon candidates can be combined into 7r°s in multiple 
ways, we select the combination with the smallest total x 2 f° r the ir° mass hypothesis. We permit 
calorimeter showers that are unassociated with charged tracks or 7r°s to remain. The summed 
energy of all showers within 90° of the lepton direction must be less than 100 MeV in each of the 
barrel and endcap region of the calorimeter. For the signal hemisphere, the energy requirement is 
increased to 200 MeV. 

The selection criteria above were designed to minimize backgrounds, which arise primarily from 
non-tau decays such as e + e~ — * qq or e + e~ — > T(45) — > BB or from feedthrough from other tau 
decay modes. A data sample which satisfies the above selection criteria on the signal side of the 
event, but with an invariant mass above the tau mass on the tag side, is used as an estimator 
of non-tau backgrounds. The background estimate is obtained by scaling the number of these 
estimator events above the tau mass on the signal side to the number of data events above the 
tau mass on the signal side minus the Monte Carlo predicted number of taus misreconstructed 
above the tau mass. The background estimates obtained in this way are consistent with those 
obtained by using e + e~ — > qq Lund |l5|] generated Monte Carlo events. Monte Carlo events with 
full detector simulation Jl6| based on the event generators KORALB llTJ and PHOTOS E8| were 
used to estimate tau decay backgrounds in the signal hemisphere (feed-across). 

The resulting event distributions, shown in Fig. |], contain 266 events in the 5-7T final state, of 
which 8 are above the tau mass. The 3ir2ir° mode contains 207 events, of which 13 are above the 
tau mass. Monte Carlo simulations show a reconstruction efficiency of (3.08 ± 0.10)% for the 5tt 
mode and (0.43 ± 0.02)% for the 37r27r mode. In both modes the number of events observed is 
consistent with the expectation from the world average branching fractions [19]. 

Each signal event is represented by a point in the two-dimensional plane formed by the hadronic 
energy scaled to the beam energy (Ex/Eb) versus hadronic mass (Mx), with the estimated errors 
represented by an ellipse. The sensitivity to neutrino mass is largest in the region near Mx = m T 
and Ex/Eb = 1, as shown in Fig. [2|. The data events are shown in Fig. ^. We determine the 
tau neutrino mass using the events lying in the fit region Q shown in Fig. ||. There are 36 5tt 
events in the fit region, with negligible tau backgrounds and an estimated 0.3 events from non-tau 
backgrounds. There are 19 3tt2tt° events in the fit region, with 1.0 of these expected as feed-across 
from other tau decay modes and 0.4 expected from non-tau backgrounds. Details of the analysis 
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TABLE I. Summary of signal, background, efficiency, resolution and upper limit by mode. 



Mode 


5tt 


3vr2vr 


Total Events 


266 


207 


Events in Fit Region 


36 


19 


Signal Region Purity (%) 


99 


93 


Selection Efficiency (%) 


3.08 


0.43 


Typical Mass Resolution (MeV/c 2 ) 


15 


25 


Typical Energy Resolution (MeV) 


25 


50 


Upper Limit @ 95 % CL (MeV/c 2 ) 


31 


33 



are summarized in Table |[ 

We extract a measurement of the tau neutrino mass using an unbinned extended maximum 
likelihood technique. The procedure fits for one parameter (m u ), taking as input the measured 
hadronic energy and mass of the events in the fit region shown in Fig. ^. The likelihood function, 
detailed in the appendix, is composed of a Poisson factor, expressing the number of events expected, 
times the sum of a signal term and a background term. The likelihood is calculated as a function 
of neutrino mass using a novel technique. Instead of using explicit parameterizations of the Monte 
Carlo in a likelihood convolution [^], Monte Carlo signal events are used in the evaluation of 
the likelihood, directly implementing the best knowledge of all physics effects including initial state 
radiation and detector acceptance. The only explicitly parameterized term is the detector smearing 
function. 

The resulting extended likelihood is shown in Fig. ||. We define^ the 95% confidence level (CL) 
upper limit by integrating defined likelihood above zero mass to its ninety-fifth percentile. We 
find 95% CL upper limits of 33, 31, and 27 MeV/c 2 for the 3tt2tt°, 5tt and combined samples. A 
separate analysis using only mass information yields a combined upper limit of 31 MeV/c 2 . 

The interpretation of this limit as a meaningful statement about probability is conditional 
upon the measured likelihood being representative of an ensemble of similar experiments. The 
event distribution and the observed number of events at the endpoint are consistent with our 
Monte Carlo estimation with zero neutrino mass, using the spectral function tuned to the data as 
described below. An ensemble of Monte Carlo experiments using statistics compatible with those 
we observe reveal that a smaller upper limit is obtained in 67% of experiments with a massless 
neutrino. The average expected upper limit given our sample size is found to be 25 MeV/c 2 . 

Following the conservative prescription used by the LEP experiments [20], a linear systematic 



8 The method for extracting an upper limit from a likelihood distribution at a given confidence 
level is not unambiguously defined; the method used here differs from that used in the analysis of 
Ref. ||. Therefore comparisons of upper limits among different experiments must be done with 
care. 
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TABLE II. Systematic error summary. 



Son vcp 


Error fMeV /r 2 ) 


Snppfyfl 1 Ki 1 n pfinn 


1 Q 

-L . £7 


Smearing Tails 


1.6 


Mass Scale 


1.1 


Smearing Parameters 


0.7 


MC Statistics 


0.6 


Background 


0.4 


Energy Scale 


0.3 


Total 


2.9 



error, A(mg5), is added to the limit. The error is defined by 

A(m 95 ) 2 = S(m 95 - m 95 ) 2 , 

where the sum is over systematic error sources, and the difference represents the change in the 
nominal upper limit from each error source. Systematic errors are summarized in Table ||. 

The largest systematic error is due to uncertainties in the spectral functions. We obtain the 
spectral functions from a fit to a statistically independent sample of decays recoiling against pionic 
decays of the tau, excluding the 100 MeV/c 2 near the tau mass and including a background term. 
Each fit is a sum of two phenomenologically inspired spectra [21-23]; the dominant spectrum is the 



distribution predicted for the 5 pion decay of the tau as derived from the CVC predictions for 4 
pion decay and soft pion theorems, and the other is the harder spectrum CVC predicts for 6 pion 
decays. The systematic variation for the spectral function is taken as the allowed variation in the 
fit parameters and results in errors of 1.5 MeV/c 2 from the 5ir mode and 1.1 MeV/c 2 from the 
3tt2tt° mode. 

The next largest systematic error comes from the modeling of the detector smearing function. 
The smearing function is approximated by a sum of three two-dimensional Gaussians to model its 
extended tails. The smearing function parameters are derived from the Monte Carlo simulation, 
while the widthes of each Gaussian are proportional to the event error, as estimated from the prop- 
agated tracking error matrix and the photon reconstruction errors. The errors from the modeling of 
the tails of the smearing function are estimated by using a fourth Gaussian component in separate 
fits to charm and B decay data. The mass (energy) smearing is studied with charm (B) decays to 
final states with similar numbers of tracks as our samples. B decays are used since the B, created 
nearly at rest in the CLEO detector frame, allows the comparison of reconstructed energy with the 
beam energy. This smearing uncertainty is 1.4 MeV/c 2 for the 3ir2ir° mode and 0.6 MeV/c 2 for 
the 5tt mode. In addition, all smearing function parameters are allowed to vary within their fitted 
uncertainties, resulting in another 0.5 MeV/c 2 uncertainty for each of the two modes. 

The uncertainty in the absolute mass and energy scales of 1 MeV/c 2 (as estimated from charm 
decays) and 2.7 MeV (from B decays), respectively, contribute 1.1 MeV/c 2 and 0.3 MeV/c 2 to the 
final systematic error estimate. 

Other sources of error considered are the Monte Carlo statistics used in the likelihood evaluation 
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and the background size. The total systematic error from all sources added in quadrature is 2.9 
MeV/c 2 , resulting in a final 95% C.L. upper limit on the mass of the tau neutrino of 30 MeV/c 2 . 

The use of the Poisson term in this analysis' likelihood is an improvement over a previous 
technique || because it reduces the variance in fitted neutrino mass, thus decreasing the probability 
of a spuriously low limit in the presence of a massive neutrino. The present analysis also uses 
measured event errors to avoid a possible bias. The smearing shape used in previous analyses H|i~0| 
depends on repeated Monte Carlo simulation of the data points, using the measured values as input 
to the Monte Carlo event generator. Such a method neglects the fact that events reconstructed 
near the endpoint stand a larger chance of being upward fluctuations from low mass events than do 
events in the middle of the accepted distribution, and thus tends to underestimate the mass and 
energy errors associated with each data point. 

In summary, using the world's largest sample of candidate r — > 5ttv t and r — > 37r27r°^ T decays, 
we obtain a 95% confidence level upper limit on m Vr of 30 MeV/c 2 . The method used properly 
takes into account the change in the number of events expected at the endpoint as a function of 
m Vr , as well as a background contribution, and a detailed detector smearing function. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent luminosity 
and running conditions. J.R. Patterson and I.P.J. Shipsey thank the NYI program of the NSF, M. 
Selen thanks the PFF program of the NSF, M. Selen and H. Yamamoto thank the OJI program of 
DOE, J.R. Patterson, K. Honscheid, M. Selen and V. Sharma thank the A. P. Sloan Foundation, 
M. Selen thanks Research Corporation, S. von Dombrowski thanks the Swiss National Science 
Foundation, and H. Schwarthoff thanks the Alexander von Humboldt Stiftung for support. This 
work was supported by the National Science Foundation, the U.S. Department of Energy, and the 
Natural Sciences and Engineering Research Council of Canada. 



APPENDIX: APPENDIX 



The likelihood is given by: 
C(m v ) - 



V (N obs ) Yl £dat{XdataWdata, m iJ, Nbgd) 
Events 



The first factor is the Poisson probability of the number of events observed in the fit region, N f, s . 
The second factor is the product likelihood over the data events observed X^ata > an d their associated 
errors (Jdata, given an expected number of background events, N^gd- 
The Poisson factor is 



V(N t 



obs) 



(N (m„) 



+ N 



bgd 



N„ 



-(N(m v )+N bgd ) 
AW 



The expected number of signal events in the fit region Q, N (m u ), is a function of neutrino mass 
and is extrapolated from the number of events seen in the control region 1Z using the partial width: 

In dT (™«- 



N(m L 



-K 



n 

data • 



J n dT (m u = 0) " 

The control region TZ is chosen away from the neutrino mass sensitive region so that the number 
of events expected in it does not depend strongly on the neutrino mass (see Fig. The partial 
width is fit from an independent data sample as described in the text. 



The event likelihood can be expanded as: 

£dat(XdataWdata, m vi Nbgd) = (1 Tjf^) C sig{X data\& data i ^v) 

obs 

+ ( at q )£bqd (Xdata\&data 

Kbs v 

The first term of this expression is the neutrino mass likelihood for a pure signal. The second term 
expresses the background shape. The likelihoods C s i g and Cbgd are normalized over the fit region 
Q and neutrino mass. The notation Xdata is to stress that we use the reconstructed (smeared) 
variables for the data. 

The background likelihood, Cbgd, is a parameterization of the background event distribution 
from the high tag-mass sample described in the text. 

The signal likelihood, £ s ig, is determined from Monte Carlo simulations, taking into account 
detector smearing, the physics of tau decays, and initial state radiation effects. The analytic 
expression for this likelihood, as used in previous analyses (see for instance ^0|), can be written as: 

r (f- I \ _ IaUX dX G(X,X D a ta,aDatg) e(X) fph ys (X\m u ) 

sig\J*- data\@ data i Tfl v ) — 



M{m u )= [ dX [ dXG(X,X,a)e(X) f Phys (X\m h 



where 

• X is the reconstructed data point coordinate: (M,E); 

• fphys(X\m u ) represents the physics, and is a convolution over the differential decay rate 
dY(m u , X)/dq 2 , and corrections due to initial state radiation, and beam energy smearing; 

• G(X,X 

; &data)i a detector smearing function, is the probability of an observed X given a true 
X and an estimated error Odata \ 

• e is the total detector efficiency 

Noting that Monte Carlo events are naturally generated with a probability density proportional to 
fphysi the signal likelihood for a given data event is rewritten as : 

,~ , Y<MC\m v =()G{X da ta-,XMC->Vdata) *VV{mv\XMc) 

J~-sig\^-data\^datai m u ) — 7T} \ • 

N(m v ) 

The sum is taken over all accepted Monte Carlo events, generated with a massless neutrino. The 
term W(m u \XMc) represents the relative weight of the Monte Carlo event for a massive neu- 
trino compared to a massless neutrino. It is calculated at the generated (M, E) rather than the 
reconstructed (M,E), and is written as: 

yv{m u \X MC ) = ^ r9 {Xmc) . 

dT [m v = 0\Xmc) 

This is simply the ratio of differential decay rates of a tau at the (generated) value of q 2 appropriate 
for this Monte Carlo event and this value of neutrino mass relative to a massless neutrino. Note 
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that the matrix element factorizes and that the hadronic currents cancel in this ratio. The function 
ensures that only Monte Carlo events which remain physical at the new neutrino mass are used, 
and is implemented by checking that the scaled energy of the Monte Carlo event falls into the 
physically allowed interval: 



[ U else 



The kinematic limits are given by: 

V± = ^2(Er/E B ) [m 2 + q 2 - ml ± (l - M 2 /El) 1 (m 2 , q 2 , u 2 )^ , 

X(x, y, z) = x 2 + y 2 + z 2 - 2(xy + yz + xz). 

The factor AT(m u ) ensures that the probability density used in the event likelihood is normalized 
for any neutrino mass: 

M{m v )= J2 ([ G(X,X MC ,a)dx)w(m u \X M c). 

The likelihood function used has been tested on Monte Carlo samples of decays with both 
massless and massive neutrinos and reproduces the input neutrino mass with no apparent bias. 
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FIG. 1. The mass distributions of the 57r (a) and 3tt2tt° (b) data samples. The dashed 
histogram represents the shape expected from Monte Carlo normalized to the number of events 
below the tau mass. The background estimate is shown as a solid histogram and is normalized for 
display purposes to five times its nominal value. 



12 



1160398-001 



LU 



1.0 
E0.9 

(0 
0) 

, m 0.8 



LU 



X0.7 



0.6 
0.5 



i i 



i i 



Control Region 








1.3 



1.4 



1.5 



1.6 



1.7 



1.8 



M (GeV / c ) 



FIG. 2. The number of events in the the fit region Q relative to the number in the control 
region 1Z in the scaled hadronic energy versus hadronic mass plane is a function of neutrino mass. 
Kinematically allowed neutrino mass contours are drawn for neutrino masses of 0, 30, 60 and 100 
MeV/c 2 . Note the typical error ellipse drawn in the fit region. 
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FIG. 3. The hadronic scaled energy versus mass distribution of the 57r (a) and 3tt2tt° (b) data 
samples in the fit region. Ellipses represent 1 a resolution errors and incorporate reconstruction 
induced systematic offsets. Kinematically allowed contours are drawn for neutrino masses of 0, 
30, 60 and 100 MeV/c 2 as solid, dashed, dot-dashed and dotted lines respectively. Events below 
the kinematically allowed region are fully consistent with signal events, once initial state radiation 
effects are considered (as they are in the likelihood.) 
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FIG. 4. The likelihood as a function of neutrino mass for the 5-zr (a), 3-7r27r° (b) and combined 
(c) data samples. The 95% CL upper limits, uncorrected for systematic errors, are at 31, 33, and 
27 MeV/c 2 , respectively. 
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